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Abstract 

Fluctuations of the color superconducting order parameter in dense quark matter 
at finite temperatures are investigated in terms of the phenomenological Ginzburg 
- Landau approach. Our estimates show that fluctuations of the di-quark gap may 
strongly affect some of thermodynamic quantities even far below and above the 
critical temperature. If the critical temperature of the di-quark phase transition 
were rather high one could expect a manifestation of fluctuations of the di-quark 
gap in the course of heavy ion collisions. 
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The quark-quark interaction in the color antitriplet channel is attractive driv- 
ing the pairing, cf. [1]. The problem has been re-investigated in a series of 
papers following Refs. [2,3], see review [4] and Refs therein. The attraction 
comes from the one-gluon exchange, or from a nonperturbative 4-point inter- 
action motivated by instantons [5] or nonperturbative gluon propagators [6]. 
The zero-temperature pairing gap A was predicted to be ~ (20 -h 200) MeV for 
the quark chemical potentials yUg ~ (350 -h 500) MeV. The critical temperature 
can be estimated as — 0.56A, as in the standard BCS theory, cf. [4]. 

One expects the di-quark condensate to dominate the physics at densities be- 
yond the deconfinement /chiral restoration transition and below the critical 
temperature. Various phases are possible. E.g., the so called 2-color super- 
conductivity (2SC) phase allows for unpaired quarks of one color. There may 
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also exist a color-flavor locked (CFL) phase [7] for not too large strange quark 
masses [8], where the color superconductivity (CSC) is complete in the sense 
that the di-quark condensation produces a gap for quarks of all three colors 
and flavors. The value of the gap is of the same order of magnitude, as that 
in the two-flavor case. 

The high-density phases of QCD at low temperatures may exist in interiors of 
most massive neutron stars affecting the cooling, rotation and magnetic field, 
cf. [9,10,11]. It is also possible to ask whether CSC is relevant for the terrestrial 
experiments? To produce a color superconducting matter in the laboratory 
one would need to cook a dense and not too hot baryon enriched matter. The 
nuclear matter prepared in heavy ion collisions at SISIOO, AGS, SPS or RHIC 
in all the cases has presumably not sufficiently high baryon density and, on 
the other hand, the matter prepared at SPS or RHIC is feasibly too hot in 
order to expect a manifestation of the CSC. The most relevant is probably the 
GSI "Compressed baryon matter" heavy ion collision facility which may cook 
enough dense and not too hot quark-gluon matter. The baryon densities up 
to ten normal nuclear matter density (< lOpo) at temperatures T < 170 MeV 
are expected to be reached at an initial collision stage. It is supposed that the 
system is in the quark-gluon plasma state at such conditions. Then the system 
expands and cools down. In this process the temperature decreases up to T ~ 
140 MeV at still rather high baryon density at an intermediate collision stage. 
Although the temperatures at relevant densities are most likely larger than 
the critical temperature Tc of CSC, one may rise the question on a possible 
manifestation of the precursor phenomena of the CSC phase transition, if Tc 
is rather high (Tc ^ (50 ^ 70) MeV). Recently, ref. [12] considered such a 
possibility within the Nambu-Jona-Lasinio model and demonstrated that the 
fluctuating pair field results in a prominent peak of the spectral function, which 
survives in the temperature interval \T — Tc\ < 0.2Tc. Besides, it is interesting 
to investigate the role of the order parameter fluctuations for T < Tc, if Tc is 
rather low (Tc ^ (30-^50) MeV), that may affect the neutrino radiation of the 
most massive hot neutron stars at an initial stage of their evolution. 

In this Letter we will study precursor phenomena of the CSC in the framework 
of the phenomenological Ginzburg - Landau approach. To be specific, we will 
consider condensates with total angular momentum J = that are antisym- 
metric in color and fiavor. Such pairing states can occur in the weak coupling 
limit because one-gluon exchange is attractive in the color anti-triplet channel. 

The Fourier component of the density of the thermodynamic potential (ther- 
modynamic potential per unit volume V) in the superconducting quark matter 
with the di-quark pairing can be written in the following form [11,13], cf. [1,10], 
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n^nn + J2i-co\drdU'' + c\Vdl,\'') + aD + ^D^ , (1) 

a,i 

The Greek indices a, /5 = {i?, 5, G} count colors, the Latin indices i = {u, d, s} 
count flavors. The expansion is presented up to the fourth order in the di-quark 
field operators (related to the gap) assuming the second order phase transi- 
tion, although at zero temperature the transition might be of the first order, 
cf. [14]. Qn is the density of the thermodynamic potential of the normal state. 
The order parameter squared is D = |rfisp = |ci_Rp + |c^Gp + MsP, c^rHc^gIMb 
for the isoscalar phase (IS), and D = 3|(icFLp, |c?Rp = |o?Gp = MsP = McflP, 
d*ji ■ da = d*Q ■ ds = d% ■ dn ^ Q for the CFL phase, da = {o^q, d'^, d^}. The 
so called 2SC phase is the IS phase with unpaired quarks of one color (to say 
R). The interaction of the di-quark field with fiuctuative gluon fields is usu- 
ally introduced in the standard way through the corresponding gauge-shifted 
derivatives. Recent paper [15] demonstrated that the gluon fiuctuations con- 
tribute essentially to thermodynamic quantities only at temperatures in a 
narrow vicinity of the critical temperature, AT/Tc ^0.1 for relevant values of 
parameters. In this small temperature interval, AT, gluon field fiuctuations 
change the nature of the phase transition (from the second to the first order). 
However the value AT is much less than the temperature region where fiuc- 
tuations of the di-quark order parameter might be important, as we will show 
below. Thereby in the latter discussion we suppress the gluon fields as well as 
the discussion of any peculiarities of this narrow temperature region near T^.. 

Near the critical point coefficients of (1) can be expanded in t = (T — Tc)/Tc. 
In the weak coupling limit they render: 



a^aoHT/Tc)^aot, ao = 6i = 62 = V(l + t) " ^^^^ 



3(1 + 7) 



, co = 3c. (2) 



For the classical (mean) fields 7 = 1 in the IS phase and 7 = 1/3 in the 
CFL phase, cf. [11], fig = /j^b/^, ^^b is the baryon chemical potential (the 
contribution of the strange quark mass is neglected), C(3) = 1.202... The 
critical temperature Tc is the same for the IS and the CFL phases for — > 0. 
For rus 0,Tc would depend on that would result in a smaller Tc for the 
CFL phase than for the IS phase. 

The value of the order parameter follows by solving the equation of motion 
for the field operators, 5(}/5dl^ — : 
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- cod'X + - ad}^ - bDdi = 0. 



(3) 



The stationary, spatially-homogeneous mean field solution of (3) (without tak- 
ing into account of fluctuations) is 



DMF = -ae{-t)/b, 50mf = -^e(-i), (4) 

where the step function 0(— t) = 1 for i < and 6(— t) = for i > 0. 

To be specific, discussing T < Tc we further consider the CFL case. For the 
finite system of a large spherical size R ^ ^, with ^ having the meaning of 
the coherence length, we obtain 



6' I 



{R-r) 




(5) 



where A^^c = 3 is the number of colors, Nf = 3 is the number of flavors and to 
be specific we assumed the simplest structure c?^ jy^p oc 5^. 

Now we will consider fluctuations. Below Tc we present o?^, as (i^ = rf^ ^ + li^', 
and above as = d^', where index "c" labels the classical solution. 

Since 5F{y,T) = 6fl{fi,T), cf. [16], the density of the free energy of the CFL 
phase expanded in (li^')^ terms renders 



5F^5Fc + 5F' = Y.i-<^o\drdlf + c\VdlJ) + aD^ + ^ilM^ 
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k k 

+ i^H^^[NcNf^mr. (6) 

k 

The last term is introduced within the self-consistent Hartrcc approximation 
(within the $ functional up to one vertex), uh = 10/9 accounts different 
coefficients in $ functional for the self-interaction terms (d^ - for the given 
field d and {d}^Y{d^pf terms), cf [17]. We presented d'^' = Efe <'^fce-^'=''^^ k^' = 

— * 

(a;, k), Xfj, — (r, — r), and introduced the notation 



E Ma fel' = E \<t>'k? = ^ I ^ ■ (7) 
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In Matsubara technique the temperature dependence is introduced by the 
replacements: uo ^ uo^ = 2ninT, -i J T En=!^oo / • We restricted 

ourselves by the self-consistent Hartree approximation. Effects of the damping 
of fluctuations, being produced by the two- and more vertex diagrams of 
are beyond the scope of this simple approximation. As we argue below, the 
main effects we discuss in this paper are not essentially affected by the width 
terms. 

Variation of eq. (6) over (f)'^ yields the spectrum of fluctuations 

coo;^ - ck^ - = 0, (8) 
with the squared mass parameter 

= r]\a\ = \{h + b2/3)D, + a + + h2)N,Nf ^ . (9) 

k 

Notice that the real physical meaning of the effective mass has the quantity 
m/ -yJcQ rather than m, as follows from (8). 

Solution of the equation for the fluctuating field presented in the coordinate 
space is characterized by the length scale I = ^^j yjf\ and by the time scale 

r = icj'^{cr\)-^l'^. Above the critical point = 0, n? ^a^Oi^y, > 0, 

and neglecting \(^'^ terms one gets the parameter ~ 1. 

Variation of (6) over (i^*^ yields the equation of motion for the classical field 



-(6i + 62/3)DeC?L,c = 0. (10) 

Only dropping the term responsible for fluctuations we flnd — Dmf, d\^^ — 
(i^ MF) ""^^ — "^MF — for T < Tc, cf. (4), (9), and we obtain ?7mf = 2/3. 

In reality fluctuations affect the classical solution and rcnormalize the critical 
temperature of the phase transition yielding TJ*'" < T^. More generally 



m'^ml^ + Sm', 5m' ^ —{h - b2/9)N,NfY,\<P',\\ (11) 

k 

for T < TJ*^". From (10), (11) one finds Dc{m' = 0) = 0, as the consequence 
of the self-consistency of our approximation scheme. As we have mentioned 
we neglected fiuctuations of the gluon fields which change the nature of the 
phase transition but yield only a small jump of the order parameter. 
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In order to demonstrate how fluctuations of the order parameter may affect 
thermodynamic quantities let us calculate the contribution of fluctuations to 

the speciflc heat density Cy — —T ^^p-^ , cf. [16]. The mean fleld contri- 



bution is 



cr = T^e(-t), (12) 

as it follows from (1), (2) and (4). Here and below we use an approximate 
expression a ~ aot valid at T near Tc, see (2). The contribution to the spe- 
cific heat from the normal quark excitations is suppressed as oc T^'^/^e^"^*^*^^/^ 
for T ^ Tc, cf. [18]. We may reproduce correct low temperature behavior 
multiplying (12) by a form-factor / ~ 2e-^^°'>/^ {TjTf/^{l - 

The fluctuation contribution F ' can be determined with the help of the func- 
tional integration 



exp{-6W) = J D(l)'exp{-6W[(j)']), (13) 
where SW is an effective potential. Using (13), (6) we obtain 

5F' = -iN^Nf J ^ln[cou;2 _ ^^2 _ ^2]_ (^4^ 

Again within the Matsubara technique one still should do the replacement uj 
ujn = 27rinT, -i J TYJ^'^oo! (|^- We dropped an infinite constant 

term in (14). However expression (14) still contains a divergent contribution. 
To remove the regular term that does not depend on the closeness to the 
critical point we find the temperature derivative of (14) (entropy per unit 
volume) : 



SS'^?^^?^N,N,^W.r (15, 

k 

The contribution of fluctuations to the specific heat density is then given by 



^ -T^iV^iV,^^. ^ . -,„.a„/r.. (16) 
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In the second line (16) we remained only the terms quadratic in fluctuation 
fields, i.e. we assumed |^^| > |^|. 

Now we are able to calculate the quantity iG {X = 0) = J2k |0a;P ^-^d its 
temperature derivative, where iG {X) is the time-ordered Green function. 
Using the Matsubara replacement a; — > = 2ninT and the relation J^niy"^ + 
n^)~^ = -cth(7ri/), or the corresponding relation between the non-equilibrium 
Green functions iG {X) and ImC*^* at finite temperature, we arrive at the 
expression 



2^0 J (271)^ 



=cth 



2T^ 



(17) 



Using that cth(?//2) = 2nB{y) + 1, where nsiy) = (e^ — 1)~^ arc Bose oc- 
cupations, and dropping the regular contribution of quantum fiuctuations we 
obtain 



d^k 



^ 2 1 f 



(18) 



_|_ 

By index "T" we indicate the thermal contribution. 



The integration is performed analytically in the limiting cases. Let m ^ T^/cq. 
Then nsiy) — e~^. In the very same approximation one has ck"^ <^ for 
typical momenta. Then 



El«-^(^)'^'-pf-;^V rn»2T^o. (19) 
^ 8m \ nc J \ T^J 

In the opposite limiting case, m <S 2T^Jcq^ there are two contributions to the 
integral (18), from the region of typical momenta cfc^ ~ vr? and from the 
region ck'^ ^ w?. They can be easily separated, if one calculates the auxiliary 
quantity — ^^f^- In the region ck"^ ~ one may use an approximation 
UBiy) ^l/y and 



^^i^[cP ^ m^] ^ -(4.c3/^)-T^. (20) 
The region ck"^ ^ yields a regular term 
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General expression for the fluctuation contribution to the specific heat is given 
by the first hne (16) and can be resolved with the help of the self-consistent 
solution of (11), (18) (or (19), (20), (21) in the limiting cases). Assuming for 
rough estimates that fluctuations can be described perturbatively and putting 
~ m|^p = r^Mplcil, from the second line of (16) we find for T <C 



7/4 



1 + 



2\t\Tr 



exp 



1(22) 



and for T > 

Vco 



87rT2|i|V2 \cj 6Te c Vc 

In the weak coupling limit forthe CFL phase t^mf = 2/3,ao/c= 67r^T^^/(7C(3)/8), 
Co = 3c. Eq. (22) holds in the low temperature limit T < 7r{\t\/3Y/^ T^, 
whereas eq. (23) is valid in the opposite limit. The first term in (23) domi- 
nates over the second one for \t\ < 0.7, i.e., in the whole region of validity of 
(23) the main term is the first one. Its singular behavior, ~ 1/^/^, is typical 
for the specific heat in the vicinity of the critical point of the second order 
phase transition, as in metallic superconductors. 

Also in the CFL phase (T < T^) there is a contribution to the specific heat of 
Goldstone-like excitations, cf. phonons in the ordinary condensed matter. For 
T ^ mp,G, where m^.o is the mass of the pseudo-Goldstone excitation, we get 



where A^^ is the number of pseudo-Goldstone modes. One can see that the 
contribution (24) is numerically small compared to those terms (cf. (22), (23)) 
we have evaluated for temperatures within the fluctuation region. 

Comparing the mean field (12) and the fluctuation (16) contributions to the 

specific heat (in the low and high temperature limiting cases one may use 
eqs. (22), (23)) we may estimate the fiuctuation temperature Tg"^ < Tc, at 
which the contribution of fiuctuations of the order parameter becomes to be 
as important as the mean field one (so called Ginzburg - Levanyuk criterion) , 
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c 



iMF 
V ' 



for T < T^. 



(25) 



Fluctuations dominate for T > T^<. For typical values Hq ~ (350 -j-500) MeV 
and for > (50 -=-70) MeV in the weak coupling limit from (25), (22) we esti- 
mate Tq^ ~ (0.6 0.8)T'c. If we took into account the suppression factor / of 
the mean field term decrease of the mass m due to the fluctuation 

contribution (cf. (11)), and the pscudo-Goldstonc contribution (24), we would 
get still smaller value of Tg^< (< O.STg). We see that fluctuations start to con- 
tribute at temperatures when one can still use approximate expressions (22), 
(19) valid in the low temperature limit. Thus the time (frequency) dependence 
of the fluctuating flelds is important in case of CSC. 

In the condensed matter physics one usually performs calculations in the high 
temperature limit. In this limit one neglects the time (frequency) dependent 

terms considering quasi-static thermal fluctuations of the order parameter. 
Then the fluctuation contribution is determined with the help of the functional 
integration 



Comparison of eqs. (27) and (14) shows that the high temperature expression 
(27) is recovered, if one drops all the terms except n = in the corresponding 
Matsubara sum over a;„ = 2'KinT in eq. (14). The contribution of n 7^ 
terms to the is suppressed in the limit CqAh'^T'^ ^ m^. Then one 

immediately arrives at the speciflc heat given by the first term of (23). Thus 
in our case one may suppress the frequency dependence of fluctuations only 
for \t\ < 1. 

Simplifying, the energy width of the fluctuation region, where the fluctuation 
effects prevail over the mean fleld ones, is usually estimated following the 
Ginzburg criterion. The probability of the fluctuation in the volume Vq is 
given hyWr. exp{-Sn{V{i)/T). It is ~ 1 for Sn{Vfi) ~ T, where 5Q(T4) (the 
contribution to the thermodynamic potential in the corresponding variables) 
is the work necessary to prepare the fluctuation within the volume Vfl. The 
minimal size of the fluctuation region characterized by an order parameter 



d ~ is ^T)^, cf. (5), (8). Thus, taking 51](Vfl) = (5fiMF(Vfi) = 



for the typical temperature Tq, when fluctuations start to dominate, we obtain 



ex.p{-SF'/T) 
5F'[d'] 




(26) 



The integration yields 




(27) 
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- 2b 3 - v^b 3"- ^^^^ 

Although the above estimate is very rough we took care of all the numeri- 
cal factors. This allows us to notice that the value estimated from (25), 
if one uses eq. (23) for Cy remaining there only the first term, is T^^ ~ 
0.5|t(rfl^<)/t(Tfl^)|^/2yG fQj, ^ ^ 3_ All the dependencies on the pa- 
rameters in expressions for T['^ and were proven to be essentially the 
same. Fluctuation region is determined by the Ginzburg-Levanyuk parameter 
Gi — |Tc — Tg'l/Tc. The larger the value Gi is the wider is the fluctuation re- 
gion. For clean conventional superconductors [19] Gi — 10~^^ 10"^^, whereas 
for superfluid He^ and in our case Gi ~ 1. 

Now we are able to discuss the validity of the self-consistent Hartree approx- 
imation in our problem. Adding the contribution to $ functional with two 
vertices produces the sun-set diagram in the di-quark self-energy. In such a 
way beyond the Hartree approximation there appears the width term in the 
self-energy (behaving as —i'ju for small u). Such a term governs the slow re- 
laxation of the order parameter in the phase transition phenomena outside 
the equilibrium. In our case (thermal equilibrium) this term can be dropped 
compared to the CqUJ^ term, which we have in the thermodynamic potential 
(1) from the very beginning, at least in both the low and high temperature 
limits. Indeed, in the low temperature limit three Green functions entering 
the sun-set diagram produce an extra exponentially small particle occupation 
factor compared to that governs the Hartree term (19). Near Tc, i.e. in the 
high temperature limit, we may remain only n — term in the Matsubara 
sum over frequencies, as we have argued, thus suppressing both the linear and 
the quadratic terms in lu. 

All above expressions for fluctuating quantities (except the vanishing of the 
pseudo-Goldstone contribution above Tc and the appearance of an extra di- 
quark decay contribution to the width) are also valid for T > Tc, if one puts 
Dc = in general expressions or m^p = a > 0, Dc = 0, r] = rjuF = 1 in the 
corresponding approximate expressions. Above we should compare the fluc- 
tuation contribution to a thermodynamic quantity with the quark and gluon 
contributions of the normal state of the quark-gluon plasma. In the weak cou- 
pling limit for the specific heat density of the quark-gluon plasma one has, cf. 
[20], 



C-^6/.?T+l^^ + ^. (29) 

The first two terms are quark contributions and the third term is the gluon 
contribution. We omitted a contribution of strange quarks which is rather 
small for T < rus and we neglected the temperature dependent effective gluon 
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mass. We also disregarded the ag corrections to the quark and ghion terms 
since such corrections were not taken into account for condensate quantities. 

Above Tc, eqs. (22) and (23) yield rather smooth functions of T except the 
region t ^ 1. The fluctuation region is rather wide since Cy ~ Cy even 
at T essentially larger then Tc. The appearance of an extra channel of the 
di-quark decay width, 27dcc^! beyond the Hartree approximation does not 
qualitatively change the situation since 7dec is a regular function of T, 7dec oc 
T — Tc, and vanishes in the critical point. Thus the applicabihty of the high 
temperature limit (one can put ujn = in the Matsubara sum) is preserved in 
a wide temperature region, T ~ Tc, [19]. Only far above Tc situation might be 
changed. The main uncertainty comes from the values of coefficients (2) which 
were derived for T rather near T^ and in the weak coupling limit. Bearing all 
this in mind we estimate the value of the fluctuation temperature T^^ ~ 2Tc 
for typical value Tc/ Hq ~ 0.1 0.3. The higher the value Tc is, the stronger is 
the contribution of fluctuations at the given ratio T/Tc. 

Thus we see that fluctuations of the di-quark gap may essentially contribute 
to the thermodynamic quantities even well below Tc and above Tc. 

Note that several simplifying assumptions were done. The coefficients (2) were 
derived for \t\ <^ 1 in the weak coupling limit neglecting the strange quark 
mass, but applied in a wider temperature region for may be not sufficiently 
large fig. Only Gaussian fluctuations were taken into account within the self- 
consistent Hartree approximation. Thus, di-quark width effects were assumed 
to be suppressed. We used the expansion of the thermodynamic potential up 
to the fourth order terms in the mean field neglecting a small jump in the order 
parameter due to gluon fluctuations. We incorporated in the thermodynamic 
(luautities only the terms which have a tendency to an irregular behavior near 
Tc, whereas above Tc the short range correlations begin to be more and more 
important with the increase of the temperature. Therefore one certainly should 
be cautious applying above rough estimates outside the region of their quan- 
titative validity. However, as we know from the experience of the condensate 
matter physics, see [19], such extrapolation equations work usually not too 
bad even for temperatures well below and above Tc. 

So far we have discussed the specific behavior of fluctuations of the order 
parameter at fixed temperature and density. There are also fiuctuations of the 
temperature and the local quark density. They are statistically independent 
quantities [16], < STSpg >= 0, and their mean squares are 



The averaging is done over the volume, P is the pressure, Va is as above 
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the volume related to the fluctuation. In the limiting case Cy -C Cy, we 
obtain < (5T)^ > /T^ ~ pq/{Nj^Cy ), where Nj^ is the number of quarks 
involved in the volume Vq. Thus, far from the critical point the contribution 
of fluctuations is suppressed as < {STy > /T^ oc 1/ \J N^- This standard 

fluctuation behavior is essentially changed for Cy > Cy . Then we get even 
larger suppression of the temperature fluctuations. For \t\ <^ 1 in the high 
temperature limit we obtain ^< {5TY > /T^ ~ |i| — > for \t\ 0, where 
we assumed that fluctuations are most probable within the typical fluctuation 
volume T4 = 47r^^/3 and ^ oc l/^J\t\ for t ^ 0. Thus, if the system is at 
the temperature T in a narrow vicinity of T^, fluctuations of the temperature 
are signiflcantly suppressed. Being formed at < — t <C 1 the condensate 
region evolves very slowly, since the heat transport is then delayed (the typical 
evolution time is roughly r oc 1/ ^J\t\) . At temperatures T outside a narrow 
vicinity of Tc, fluctuations of the temperature resulting in a signiflcant decrease 
of the temperature in the volume T4i are rather probable, ^ < > /T^ ~ 1. 

This is the consequence of a very short coherence length and, thus, not too 
large number of quarks contained in the volume Va, ^ ^ 0.13/(Tcy^) being 
~ 0.5 fm for ~ 50 MeV and for \t\ ~ 1. Fluctuations of the quark density 
are also large for typical \t\ ~ 1, due to the smallness of Va. Thereby we 
argue that the system may produce the di-quark condensate regions of the 
typical size ^, thus feeling the possibility of the phase transition even if its 
temperature and density are in average rather far from the critical values. 

Anomalous behavior of fluctuations might manifest itself in the event-by-event 
analysis of the heavy ion collision data. In small {I) size systems, I < ^, (zero 
dimension case would be / <^ the contribution of fluctuations of the order 
parameter to the speciflc heat is still increased, see [19]. The anomalous behav- 
ior of the speciflc heat may affect the heat transport. If thermalization hap- 
pened at an initial heavy ion collision stage, and a large size system expands 
rather slowly, its evolution is governed by the approximately constant value 
of the entropy. Due to a large contribution to the speciflc heat (and to the en- 
tropy) of the di-quark fluctuations, an extra decrease of the temperature may 
occur resulting in an essential slowing of the flreball expansion process (due to 
smaller pressure). Having the di-quark quantum numbers, fluctuations of the 
gap may affect the di-lepton production rate from the quark-gluon plasma, as 
it has been noted in [12] . 

We would like once more to emphasize that the coefficients of the thermody- 
namic potential (1) obtained in the weak coupling limit might be essentially 
modified, in case if we applied the results to the fireball produced in a heavy 
ion collision, somehow changing our conclusions. E.g., the quark chemical po- 
tential decreases with the temperature. It results in an additional increase of 
the contribution of fluctuations, cf. dependences of the coefficients (2) on Hq. 
Thus we may need to analyze the strong coupling limit instead of the weak 
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coupling limit we discussed. If Hq becomes smaller than the strange quark 
mass we come from the possible 3SC phases to the 2SC phases. A more gen- 
eral discussion of modifications beyond the framework of the weak coupling 
limit can be found in [21]. A discussion how the coefficients (1) may in general 
vary within the Ginzburg - Landau approach is given in [11]. Our above argu- 
ments are admittedly speculative, and mean only to demonstrate a possibility 
of their application to heavy ion collisions and to an initial stage evolution of 
the neutron star. A more detailed treatment of the problem definitely needs a 
further work. 

Other quantities associated with second derivatives of the thermodynamic po- 
tential are also enhanced near the critical point demonstrating typical 1 / ^J\t\ 
behavior, cf. [18]. However numerical coefficients depend strongly on what 
quantity is studied. E.g. fluctuation contributions above Tc to the color dia- 
magnetic susceptibilities 



Xa^-(d'SF/dHl) Ha = cur\Aa, (31) 

are proven to be -C 1 everywhere except very narrow vicinity of the critical 
point. In spite of a smallness, as we know, in ordinary clean superconduc- 
tors the fluctuation diamagnetism turns out to be of the order of the Pauli 
paramagnetism even far from the transition [19]. 

Concluding, within the Ginzburg - Landau approach we estimated the fluctu- 
ation energy region at the GSG phase transition. The quantitative estimates 
are based on the values of parameters derived in the weak coupling limit for 
fluctuations of the GFL order parameter. Qualitative results survive also for 
fluctuations of other possible phases. We found that the frequency dependence 
of fluctuations is important for CSC in a wide temperature region. Fluctua- 
tions may contribute essentially to the spcciflc heat even at T rather far below 
and above Tc. We estimated Tfl^< < 0.5Tc and Tfl^> ~ 2Tc. Our rough estimates 
show that the high temperature GSG could manifest itself through fluctuations 
of di-quark gap in the course of the heavy ion collisions at S1S200, if the critical 
temperature of the phase transition is indeed rather high, Tc ^ (50^70) MeV. 
However quantitative results depend on the values of several not sufficiently 
known parameters and further studies are still needed to arrive at definite 
conclusions. 
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